Substituting the main group element in cobalt - iron 
based Heusler alloys: Co2FeAli x Si x . 

Gerhard H. Fecher and Claudia Felser 

Institut fur Anorganische Chemie und Analytische Chcmic, 
Johannes Gutenberg-Universitat Mainz, D-55099 Mainz, Germany 

E-mail: fecherOuni-mainz.de 
Abstract. 

This work reports about electronic structure calculations for the Heusler 
compound Co2FeAli_ a; Si a ;. Particular emphasis was put on the role of the 
main group clement in this compound. The substitution of Al by Si leads to 
an increase of the number of valence electrons with increasing Si content and 
may be seen as electron-doping. Self-consistent electronic structure calculations 
were performed to investigate the consequences of the electron doping for the 
magnetic properties. The series Co2FeAli_ a; Si :c is found to exhibit half-metallic 
ferromagnetism and the magnetic moment follows the Slater-Pauling rule. It is 
shown that the electron-doping stabilises the gap in the minority states for x = 0.5. 
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1. Introduction 

Heusler compounds [1-3] have attracted scientific and technological interest for their 
potential use as materials for magneto-electronic devices. Reason is the exceptional 
electronic structure found in many of those compounds, in particular in those based 
on cobalt. They exhibit a complete spin polarisation at the Fermi energy (e^), that 
means they behave like a metal for electrons of one spin direction and like an insulator 
for the other. Kiibler et al [4] recognised that the minority-spin density at the 
Fermi energy nearly vanishes in the Heusler compounds Co2MnAl and Co2MnSn. 
The authors concluded that this should lead to peculiar transport properties in these 
compounds because only the majority density contributes. Materials with a complete 
spin polarisation at ep are called half- metallic ferromagnets [5] , even though there do 
exist more complicated cases as classified in Ref. [6]. 

The Heusler compounds are usually ternary 2-1-1 compounds. They consist for 
the most part of two transition metals (X2, Y) and one main group (Z) element 
crystallising in the L2± structure (space group F m3m). Besides ternary X2YZ 
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compounds, there exist also large assortments of substitutional quaternary alloys of 
the type X 2 Y 1 _ X Y' X Z or X 2 YZ 1 _ X Z , X . 

One of the early substitutional series that attracted interest as potential material 
for magneto-electronics was Co 2 Cri_ a; Fe 2: Al [7-10]. Drawback of this series is that it 
is hard to be stabilised in the L2\ structure. Mostly a mixture of atoms in Y and Z 
positions is observed leading to £?2-like disorder [11]. However, the disorder destroys 
the half-metallic properties [10]. Recently, the series of Heusler alloys Co2Mni_ x Fe x Si 
has attracted particular interest because it exhibits the L2\ order over the whole 
range of x [12]. The Curie temperatures of the end members are 985 K [13,14] 
and 1100 K [15, 16] for the Mn and Fe containing compounds, respectively. The 
end members of the series C^Mni-^Fe^Si, that are the purely Mn or Fe containing 
compounds, have been used for fabrication of magnetic tunnel junctions [17, 18]. 
The tunnel magneto- resistance (TMR) ratios of 159% in the Mn compound at low 
temperature and 41% in the Fe compound at room temperature suggest that still an 
improvement in the materials is necessary for successful use in devices, in particular 
with respect to their temperature behaviour. Recently, Tekuda et al [19] reported 
about tunnel junctions build from the iso-electronic compound Co2FeAlo.5Sio.5. The 
junctions exhibited TMR ratios of 76% at 300 K and 106% at 5 K for the B2 structure 
while that with L2\ structure showed 51% and 78% at 300 K and 5 K, respectively %. 
These values of the TMR ratio are larger than the ones found using pure Co2FeAl or 
Co2FeSi electrodes. 

The temperature stability of the minority gap is one of the main challenging 
questions for materials to be used in applications. From the viewpoint of the 
electronic structure, several different effects may destroy the half-metallicity at finite 
temperatures, depending on the situation of ep. At T > 0, quasi-particle states - 
occurring close to the minority band edges - may be induced in the gap [20], for 
example by magnon excitation. In particular, a spin-rotation [21] may destroy or at 
least reduce the size of the gap. This has the effect that ep - being initially situated 
at the top of the valence band - does not fall any longer inside of the gap at elevated 
temperature. If ep is located at the bottom of the conduction band then the half- 
metallicity will be immediately lost for T > due to the occupation of minority states 
through the Fcrmi-Dirac distribution. (Note: this effect cannot appear the same way 
for ep being situated at the top of the valence band, as there are no states to be 
occupied thermally inside of the gap but only above.) Finally, the lattice parameter 
and defect densities will be changed at elevated temperatures resulting also in changes 
of the electronic structure. An already small increase of the lattice parameter may be 
able to destroy the half-metallicity if the Fermi energy was initially at T = located 
close to one of the edges of the minority gap. At the same situation, a smearing of 
the states close to ep by an increase of the defects with temperature above K, will 
also destroy the half-metallic character. Taking all those facts together, it is expected 
that a location of cf close to the middle of the minority gap will result in the most 
robust half-metallicity, supposed the gap is not too small (« 1 eV). 

Low magnetic moment compounds like Co2CrAl exhibit a variety of majority ri- 
bands crossing the Fermi energy in rather all directions of the Brillouin zone. The high 
magnetic moment C02YZ compounds (m > 4/ig) exhibit only few, strongly dispersing 
majority d-bands crossing ep mainly along high symmetry directions. These few bands 

X The TMR ratio is 175% at 300 K for optimised junctions with L2\ structure, private communication 
by K. Inomata (Tsukuba, Japan). 
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may be in favour for coherent tunnelling [22]. Finally, a mixture of the 3d elements 
on the Y position will cause different localised moments on different sites. This might 
lead to instabilities of the half-metallic character. All those aspects, fixing the rf-state 
element on the Y position and thus the kind of localised moment, location of the Fermi 
energy with respect to the gap in the minority states, and simplicities of the majority 
bands crossing it, were of prime importance for selection of the C^FeAli^Siz series 
for the present study. 

2. Computational details 

The electronic structure of the series of alloys was calculated by means of the full 
potential linearised augmented plane wave (FLAPW) method. For Co 2 FeAl 1 _ a ,Si a; , 
the calculations were carried out using the FLAPW method as implemented in 
WiEN2k provided by Blaha et al [23,24]. The exchange-correlation functional was 
taken within the generalised gradient approximation (GGA) in the parameterisation 
of Perdew et al [25]. In addition, the LDA+/7 method [26] was used to respect on-site 
correlation at the 3d transition metals, ft should be mentioned that the +U was used 
here together with the GGA rather than the LSDA parameterisation of the exchange- 
correlation functional. However, no significant differences were observed using either 
of these paramctcrisations. In WiEN2k, the effective Coulomb-exchange parameter 
U e fj = U — J is used, where U is the Coulomb part and J is the exchange part. The 
use of U e fj suppresses multipole effects. That means, it neglects the non-spherical 
terms in the expansion of the Coulomb interaction. In particular, the values for U e ff 
were set to Uc — 0-14 Ry, and Up e = 0.132 Ry, independent of the Si concentration. 
These values are able to explain the the magnetic moment in Co2Mni_ 2: Fe a ;Si over 
the whole range of Fe concentration, as was found in previous calculations [12]. Uc 
and U Fe are close to the values for the Coulomb interaction Udd for d electrons in the 
elemental 3d transition metals reported in Ref. [27]. Finally, a 25 x 25 x 25 point mesh 
was used as base for the integration in the cubic systems resulting in 455 fc-points in 
the irreducible wedge of the Brillouin zone. No noticeable changes in the precision 
of the magnetic moments or in the position of the Fermi energy were observed if 
comparing to a smaller 20 x 20 x 20 mesh. The energy convergence criterion was set 
to 10~ 5 Ry and simultaneously the criterion for charge convergence to 10~ 3 electrons. 
This combination resulted in final values being about one order of magnitude lower 
for both criteria. 

The properties of the pure Al or Si containing compounds were calculated in 
F m3m symmetry using the lattice parameter [o,ai — 5.706 A, a$i — 5.633 A) as 
found from a structural optimisation. The values are close to the experimental ones: 
oai = 5.727 Aand a,si — 5.64 A. The larger deviation for the Al compound may 
be caused by the fact that this compound does frequently not have an ordered L2\ 
structure in experiments. Following Vegard's law, a linear variation of a was assumed 
for the mixed compounds. All muffin tin radii were set to nearly touching spheres. The 
mixed compounds with x = 1/4 and 3/4 were calculated in P m3m and for x = 1/2 
in P A/mmm symmetry, similar as for the mixed Cr-Fe compounds as reported in [9]. 

3. Results and Discussion 

In the following the electronic and magnetic structure of the series Co2FeAli_ :c Sl c will 
be discussed. The end members Co2FeAl and Co2FeSi were already presented in detail 
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Figure 1. Spin resolved band structure of Co2FeAl. 

Compared are the band structures calculated in the LSDA-GGA (a,b) and the 
LDA+t/ (c,d) approaches. 

in previous work. In Rcfs. [15, 28] it was shown, that it is necessary to include the 
on-site correlation in the calculations for Co2FeSi in order to explain the experimental 
data and to find the half-metallic ground state. Details of the band structure for that 
system are found in Refs. [12,28]. Already in pure LSDA-GGA calculations, Co2FeAl 
became a half-metallic ferromagnet [9]. It can be expected, however, that on-site 
correlation plays also an important role in the Al containing compound if it does in 
the case of Si. Therefore, the electronic structure of Co2FeAl was recalculated using 
the LDA+£/ method as described in section 2. Figure 1 compares the spin resolved 
band structure of Co2FeAl calculated in the LSDA-GGA and the LDA+/7 approach. 

From Figure 1, it is seen that the inclusion of U e ff in the calculation does not 
cause pronounced changes of the majority bands. Even the flat band at about -4 eV 
below the Fermi energy is shifted by only 200 mcV to higher binding energies. This 
is remarkable as this band is mainly responsible for the localised moment at the Fe 
atom. The major impact of the Coulomb parameter is on the minority bands and in 
particular on their unoccupied part. The gap is clearly opened up and the flat, lowest 
conduction bands at the T-point are shifted by about 1 eV to higher energies. 

Figure 2 compares the spin resolved density of states for the complete series 
Co2FeAli_ a; Si a ; The calculations were carried out using LDA+(7. For x = the low 
lying s-states are found at energies below -6 eV. With increasing Si content a new 
group of s-states appears that is found for x = 1 at energies below -9 eV. In the range 
of intermediate Si concentration both groups appear. Those low lying s-states are 
separated from the p and ci-states by the Heusler-typical hybridisation gap. This gap 
is considerably larger in the Si compound (« 1.5 eV) compared to the Al compound 
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(« 0.5 cV) indicating the stronger hybridisation. This stronger hybridisation makes 
the Si rich alloys more stable compared to the Al rich part. The p-states are in all cases 
found at the bottom of the high lying valence bands above that gap. With increasing 
x, both - the majority as well as the minority channel - do not exhibit pronounced 
changes of the d-state derived densities. Their general shape stays rather unaffected. 
However, the (i-band width increases from about 5.2 eV to 6.7 eV with increasing Si 
content. At the same time, the high majority density of the localised rf-states found 
at -4 eV in the Al compound shifts to -5 eV in the Si compound. Fixing the Fermi 
energy in the minority gap at the position where it is found in the Al compound would 
thus lead to a rather unphysical enlargement of the exchange splitting. However, the 
shift of the majority density is compensated by a shift of the minority density with 
increasing Si content to lower energies. As a result, one observes a virtual movement 
of ep through the gap in the minority states. Throughout the whole series from x = 
to 1 , the band gap about ep is clearly revealed in the minority density. The properties 
of the gap are further discussed after a brief discussion of the magnetic moments. 
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Figure 2. Spin resolved density of states of Co2FcAli_ I Si I . 

The panels (a, ... , e) show - from top to button - the DOS with increasing amount 
of Si for x = 0,0.25,0.5,0.75, and 1. The DOS is calculated using LDA+l/. 



The magnetic properties are compared in Tabel 1. The nearly half-metallic 
state found for Co2FeAl already in the LSDA-GGA calculation results in a nearly 
integer spin magnetic moment (deviation smaller than 10 -3 ). This reflects the well- 
expected fact that a magnetic moment being compatible to the Slater-Pauling rule 
(m = (N v — 24) ii B , with N v being the number of valence electrons in the primitive 
cell containing four atoms [29,30]) must not result unambiguously in a half-metallic 
state. The substitution of Al by Si results in a pronounced deviation from the Slater- 
Pauling-like behaviour, in the LSDA-GGA approach. This is clear from the shift of 
the minority band-gap away from the Fermi energy in the LSDA-GGA calculations 
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(see Fig. 3(a)). 

Table 1. Total magnetic moments of ordered C^FeAli—^SL,;. 
All moments were calculated for the given super-cells. Their values are in hb and 
respect 4 atoms in the cell for easier comparison. 
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Co 4 Fe 2 AlSi 


1/2 


5.37 
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Co 8 Fe 4 AlSi 3 


3/4 


5.49 


5.75 


Co 2 FeSi 


1 


5.53 


6.0 



Other than in pure LSDA-GGA, the LDA+J7 calculations reveal clearly a linear 
dependence of the spin magnetic moment m(x) on the Si concentration x and thus on 
the number of valence electrons. The small deviation from a half-metallic state for the 
x = and 1 compounds seen in Fig. 3(b) does not lead to a discernible deviation from 
the expected integer moments of 5 /j,b or 6 /j,b, respectively. This only small deviation 
is caused by the fact that only very few states contribute to the minority density close 
to ep (see Fig. 2). Finally, it has to be noted that the LSDA-GGA calculations do 
not reflect the experimentally found magnetic moment of 6 hb f° r Co 2 FeSi, either at 
the optimised or at the experimental lattice parameter. It was carefully checked that 
this effect is not caused by a missing orbital magnetic moment. Using LSDA together 
with spin-orbit coupling and Brook's orbital polarisation resulted in m s = 5.536/is 
and mi = 0.189yU_B. This means that m to t — 5.725/i_b is still far below the value 
determined by magnetometry (this magnitude of the total magnetic moment m to t was 
observed using WiEN2k as well as the relativistic Munich-SPRKKR, not reported 
here.). Only the LDA+C/ scheme was able to reflect the value in the correct order, 
independent whether spin-orbit interaction was respected or not. 

Figure 3 compares the behaviour of the gap in the minority states of 
Co 2 FeAli_ 2: Si a ;. Shown are the extremal energies of the states involving the minority 
band gap that are the accompanied valence band maximum and conduction band 
minimum. In the LSDA-GGA approach, the small gap in the minority states is 
moved away from the Fermi-energy with increasing Si content and the half-metallicity 
becomes destroyed. Using LDA+J7, the gap has a nearly constant width of 760 meV 
over the complete series from x = to 1. From Figure 3, it is seen that the end- 
members are just at the borderline to half-metallic ferromagnetism. Starting from 
x = 0, the Fermi energy moves from the top of the valence band to the bottom of 
the conduction band at x = 1. For Co 2 FeAlo.5Sio.5, the Fermi energy is located close 
to the middle of the band gap in the minority states. It should be noted that the 
values for U e ff used here are the borderline cases for the half-metallic ferromagnetism 
over the complete series Co 2 FeAli_ a; Si a ;. Small variations of the Coulomb parameter 
U e ff will change the behaviour of the density of states as reported in Refs. [15,28] 
for Co 2 FeSi. About 10% higher or lower values will make one of the end members a 
clear half-metal and destroy at the same time the half-metallic character of the other 
one. From the theoretical point of view, the compounds Co 2 FeAl and Co 2 FeSi are 
thus very unstable half-metallic ferromagnets, if at all. 



Substituting the main group element in cobalt - iron based Heusler alloys: Co 2 FeAli^ x Si x .7 



0.8 - _ _o - cbm - - 




0.00 0.25 0.50 0.75 1.00 0.00 0.25 0.50 0.75 1.00 
Si concentration x 



Figure 3. The gap in the minority states of Co2FeAli_, c Si x . 
Compared arc the positions of the valence band maximum (VBM) and the 
conduction band minimum (CBM) as calculated by means of LSDA-GGA in (a) 
and LDA+l/ in (b). 

4. Summary and Conclusions 

The electronic structure of the substitutional series of the quaternary Heusler 
compounds Co2FeAli_ x Sia; was investigated by means of band structure calculations 
using the LDA and LDA+U approximations. It was found that the Co2FeAli_ :c Sl c 
series of compounds exhibits half-metallic ferromagnetism if using the LDA+C 
scheme. Moderate Coulomb-interaction parameters of less than 2 eV were used. For 
the two end-members, Co2FeAl and Co 2 FcSi, the Fermi energy is close to the band 
edges of the minority states. The high densities at those band edges make the half- 
metallic character of both compounds rather unstable at finite temperatures above 
K. This might be one reason explaining the low tunnelling magneto resistance ratio 
found in those compounds at room temperature. For x w 0.5, the calculations predict 
that the Fermi energy is located in the middle of the gap of the minority states. 
This behaviour will make Co2FeAlo.sSio.5 stable against temperature variations as 
discussed in the introduction. Experiments were started to verify over what range of 
compositions the series Co2FeAli_ a; Si a ; crystallises in the required L2\ structure and 
to find the most stable half-metallic ferromagnet in this series. 

In summary it was shown that the variation of the main group element in Heusler 
compounds is a strong tool in order to tune their physical properties. 
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